############
# Get data #
############

# z_star grid & transition matrix 
z_star_grid  = DataFrame(CSV.File(joinpath(ddat,"z_star_grid.csv"); header=1))[:,1]
size_z = length(z_star_grid)
transition_z    = Array(DataFrame(CSV.File(joinpath(ddat,"transition_matrix_z.csv"); header=1)))

# baseline aggregates 
baseline_aggregates_df = DataFrame(CSV.File(joinpath(ddat,"baseline_aggregates.csv"); header=1))

# load baseline parameters 
baseline_parameters_df = DataFrame(CSV.File(joinpath(ddat,"baseline_parameters.csv"); header=1))
fixed_cost_param_df = DataFrame(CSV.File(joinpath(ddat,"fixed_cost_param.csv"); header=1))

# load profits_grid_baseline
expected_profits_baseline = DataFrame(CSV.File(joinpath(ddat,"profits_grid_baseline.csv"); header=1)).profits_expected
profits_NC_baseline = DataFrame(CSV.File(joinpath(ddat,"profits_grid_baseline.csv"); header=1)).profits_NC

# get baseline unconditional z_star distribution 
z_star_uncond_distrib = DataFrame(CSV.File(joinpath(ddat,"z_star_uncond_distrib.csv"); header=1)).share

# get baseline stationary distribution 
baseline_SS_distrib = DataFrame(CSV.File(joinpath(ddat,"baseline_SS_distribution.csv"); header=1)).share
